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ABSTRACT 

We present the results of a high-resolution spectral differential imaging sur- 
vey of 12 nearby, relatively young field L dwarfs (< 1 Gyr) carried out with 
HST/NICMOS to search for planetary mass companions at small physical sep- 
arations from their host. The survey resolved two brown dwarf binaries: the 
L dwarf system Kelu-lAB and the newly discovered L/T transition system 
2MASS 031059+164815 AB. For both systems common proper motion has al- 
ready been confirmed in follow-up observations which have been published else- 
where. The derived separations of the binaries are smaller than 6 AU and con- 
sistent with previous brown dwarf binary statistics. Their mass ratios of q > 0.8 
confirm the preference for equal mass systems similar to a large number of other 
surveys. 



1 This work is based on observations made with the NASA/ESA Hubble Space Telescope, obtained at the 
Space Telescope Science Institute (STScI) and associated with program GO-10208. STScI is operated by 
the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-26555. 
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Furthermore, we found tentative evidence for a companion to the L4 dwarf 
2MASSW 033703-175807, straddling the brown dwarf/planetary mass boundary 
and revealing an uncommonly low mass ratio system (q w 0.2) compared to the 
vast majority of previously found brown dwarf binaries. With a derived minimum 
mass of 10 Mj up to 15 Mj up a planetary nature of the secondary cannot be ruled 
out yet. However, it seems more likely to be a very low mass brown dwarf 
secondary at the border of the spectral T/Y transition regime, primarily due to 
its similarities to recently found very cool T dwarfs. This would make it one 
of the closest resolved brown dwarf binaries (0.087" ± 0.015", corresponding to 
2.52 ± 0.44 AU at a distance of 29 pc) with the coolest (T cff « 600-630 K) and 
least massive companion to any L or T dwarf. 

Subject headings: planetary systems — stars: low-mass, brown dwarfs — bina- 
ries: individual (2MASSW 033703-175807) — techniques: high angular resolution 



1. Introduction 



Since the discovery of the first extrasolar planet around 51 Pegasi by iMayor fc Queloz 
in 1995, more than 480 extrasolar planets (exoplanets) in more than 410 planetary systems 
have been identified as of August 20100. However, most of these detections were made 
using indirect detection methods like radial velocity (RV) surveys, the transit method or mi- 
crolensing events. Further, the RV technique can in general only provide minimum masses 
(M sin i) and thus limited knowledge of the astrophysical properties of the companion. From 
the beginning, these exoplanets revealed a remarkable diversity in their physical and orbital 
characteristics like separation, eccentricities and masses and thus challenged the planet for- 
mation scenario. In particular, the period-mass distribution shows tendencies which are hard 
to explain within the framework of current models (e. g. massive "hot Jupiter" were found at 
very small orbits, even down to a semi-major axis of only 0.02 AU). A revision of the stan- 
dard scenario seemed to be required. However, the discovery space of R V searches does no t 



yet cover planets at distances larger than about 5AU (for statistics see lMarcy et al.ll2003l ). 
The question is, if the observed period-mass distribution is, at least partly, a selection effect 
(higher sensitivity at shorter orbital periods) and thus might provide only an incomplete pic- 
ture of the overall planet population, or if giant planets are really rare at larger separations. 
In addition, the indirect RV and microlensing techniques do not enable any photometric or 



1 http : / /exoplanet . eu/cat alog . php 
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spectroscopic information of the exoplanet companions and their physical properties such as 
brightness, color, effective temperature or composition. 

Therefore, a lot of attention has been directed towards direct imaging of planetary 
mass objects (PMO) in the last years. A direct detection would enable extensive follow- 
up studies to characterize the planets, but the main challenge of this method is to deal 
with the huge brightness contrast and the small angular separation between the faint planet 
companion and its host star. With the possibility of high-resolution space-based observations 
and the improvement of adaptive optics (AO) systems on large ground-based telescopes, the 
necessary sub-arcsecond spatial resolution became achievable. In addition, planetary mass 
objects and brown dwarfs are hotter at younger ages and hence the brightness contrast 
between the primary and a possible companion will be much smaller. Thus, in order to 
push the direct detection threshold down to the planetary mass regime, systematic searches 
focused especially on young stars in nearby associations. 

First extensive Hubbl e Space Telescope (HS T) programs to search for substellar compan- 
ions to post-TTauri stars (IBrandner et al.ll2000l). as well as to v ery low mass (VLM) stars and 
brown dwarfs in the Pleiades open cluster (IMartm et al.ll2000l ) did not reveal any companion 
at intermediate orbits. Since then, several other extended HST surveys of free-floating ultra- 
cool dwarfs and brown dwarfs showed that roughly 20% of the m have substellar companions. 



Almost all of them constitut e close to equal-mass systems (e. g. lBouy et al.ll2003l ; iReid et al. 
200ll ; iBurgasser et al.ll2003bl ) and hence, rather resemble more massive, stellar binaries than 
planetary systems. Consequently, some interesting questions arose: 

Do "single" brown dwarfs evolve similar to single stars? 
Are they capable of forming planetary systems? 

The detection of circumstellar disks around youn g brown dwarfs revealed that they at 
least posses the raw material to form pl anets (e. g.[. 



Liu et al. 



2003 



Luhman et al1 l2006aL I2OO8I : IScholz et al" 



Apai et al.l 120021 ; iPascucci et al 



Javawardhana et al.ll2003l; later confirmed by result s of e. g.lApai et al 



2008, see also reviews in Luhman et al. 



2003; 
2005; 



1200 



3 and 



Hennindl2008f l. 

On that account, we started in 2004 the first spectral differential imaging (SDI) survey with 
HST to directly detect planetary mass objects in close orbits around nearby free-floating 
brown dwarfs. A direct image detection would have important implications on our under- 
standing of planetary mass objects, their origin, the formation of substellar bodies in general 
and binarity in the brown dwarf domain in particular. 



Brown dwarf multiplicity properties exhibit significant differences in frequency, sepa- 
ration and mass ratio distribution compared to Sun-like stars. While the binary fraction 
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for Sun-like stars is ~57% (jDuquennoy fc Mayor! Il99lf). this fraction steadily decreases to 
« 26-42% for early- to mid-type M dwarfs (e.g. iReid &: Gizid 119971 ; iBergfors et 
down to ~10-30% for the very low mass stars and brown dwarfs (e.g. iBouy et al. 



2003 



Burgasserl 120071 ; iJoergend 120081 ; I Goldman et al.l 120081 ) . The same steady decline is true for 
the semi-major axis distribution w hich is spread over a wide range of separations for G 
dwarfs with a peak around 30 AU (jDuquennoy fc Mayor! Il991l). smaller separation s for M 
dwarfs where the distribution peaks between 4 and 30 AU ([Fischer fc Marcylll992l ). and a 
very narrow separation distribution for the brown dwarfs which peaks at 3- 10 AU (e.g. lAllen 



20071 ). The mass ratio distribution for very low mass stars and brown dwarfs reveals a clear 
tendency to equal mass system components with a distinct peak around In contrast, 

low q values a re much more common at all separations for binary systems of solar-type stars 
(see review by lBurgasser et al.l 120071 and references therein). 

Recent analysis of radial velocit y surveys confirm a 6- 8 % fr eque ncy of planetary mass 



companions to Sun-like stars (e.g. iGrether fc Lineweaverl 120061 and iMordasini et al. 



2009 



with references therein) at relatively small separations which peak around 1 AU (jUdry fc Santos 



20071 ). While commonly observed mass ratios between solar type stars and substellar (plan- 
etary) mass objects range between 1:1000- 1:100 (cf. Extrasolar Planets Encyclopedia^) one 
would expect a much higher mass ratio of 1:20-1:5 for a system consisting of a planetary 



mass object orbiting a brown dwarf due to the lower initial mass of the primary. lApai et al. 



( 120081 ) point out that depending on the planet formation process, the frequency for lower- 
mass planets is expected to be higher than for solar type stars. For the semi-major axis 
distribution of such a system one would also expect relatively small separations of a few 
AU due to the much smaller sizes of the circumstellar disks around young brown dwarfs 
compared to solar type stars. 



2. Observing strategy 



Because of the higher achievable contrast and the lack of suitable AO reference stars, 
the observations were conducted from space with HST as explained in the following. The 
observing strategy that suited best to reach our goal of detecting planetary mass objects 
around brown dwarfs and to overcom e the remain i ng br ightness contrast, was the spectral 
differential imaging (SDI) technique (jSmith et al.lll987l ). This technique takes advantage 
of the fact that a cool (< 1300 K) giant planet or low-mass brown dwarf (i.e.T dwarf) is 
much fainter in certain molecular bands due to strong absorption than in the neighboring 



2 lit t p : / /exoplanet . eu / catalog . php 



-5 - 



continuum. Already lAllard fc Hauschildtl (119951 ) stated that the spectral energy distribution 
(SED) of brown dwarfs is very peculiar. The molecular opacities, which globally define the 
continuum of brown dwarfs, cause the SED to peak around 1.1 /im for solar metallicities, 
almost independent of their effective temperature. Further, the observed J-K colors for L 
dwarfs are red (about +2 mag), while they are blue (-0. 1 to +0.4 mag) for the T dwarfs. 
This is caused by dust grain formation ( lAllard et al.ll200ll ). Therefore, the brightness differ- 
ence of any L dwarf - T dwarf pair is smallest around 1.1 /mi, making the J -band the best 
spectral regime to detect T - type or even cooler companions to L dwarfs. 

Figure [1] shows the spectral features of L and T dwarfs in this wavelength regime in 
detail. It reveals that L dwarfs have a comparative flat SED between 1.06 /mi and 1.15 //m, 
while the flux of a T dwarf drops significantly at the beginning of the molecular water 
absorption band at ~ 1.10 mn. The T dwarf p lot in this Figure [T] additionally includes 
a synthetic spectrum from iBurrows et al.l (120031 ) for a 5Mj up , 1 Gyr old planetary mass 
object with T e g = 312 K and shows, that this effect is still present at these very low effective 
temperatures. Thus, a subtraction of two images obtained simultaneously, one in and one 
off, yet near this molecular band, will cancel out most of the equally bright PSF structure 
of the primary (L dw arf), revealing the much fainter signal from a cool, possibly planetary 
mass companion (see iRosenthal et al.l Il996l ; iRacine et al.lll999l ). The narrow-band filters 



F108N and F113N of the NICMOS1 camera (NIC1) on board HST match particularly well 
the relevant bands (as shown in Figure [1]) and are ideally suited for the differential imaging 
technique. 



3. Sample selection 

Given the angular separation and luminosity ratio problem described above, the list of 
targets was assembled with the focus on brown dwarfs in proximity to the Sun and with a 
relatively young age. A closer distance to the Sun translates in a larger angular separation for 
a given physical separation between the brown dwarf and an eventual planetary-mass object, 
thus companions on smaller orbits can be detected. In addition, the detection threshold 
caused by sky background decreases for a fainter planetary mass object, since its apparent 
brightness becomes larger, the closer it is to the Sun. 

Concerning the age, priority was given to targets with estimated ages younger than 1 
Gyr. Giant planets, as well as young brown dwarfs are more luminous when they are young 
due to the remaining gravitational contraction energy from their formation process. Left 
without a resource of producing new energy they cool down rapidly and become dimmer with 
time. Based on evolutionary models and synthetic spectra, a young (100 Myr) giant planet or 
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brown dwarf is ~ 100 - 20 times more self-luminous than a 5 Gyr old one (e. g. lBurrows et al. 



19971 ; iBaraffe et al.ll2003l ). Hence any planetary mass companion to an old brown dwarf would 
be relatively cool and too faint for detection. 

The age estimation of the targets in this survey was not straightforward since all of them 
are field br own dwarfs withou t any hints of being a member of an association or moving group. 



In general, I Allen et al.l (120051 ) showed that proceeding down the spectral sequence from L to 
T dwarfs, the average age decreases and the relative proportion of young dwarfs (<2 Gyr) 
increases if compared to M dwarfs. This is mainly caused by the lack of stable hydrogen 
burning in the bro wn dwarf core and th erefore their rapid cooling through these temperature 



regimes with time (IBurrows et al.ll200ll ). Our targets were selected according to spectral age 
indicators, such as the presence of lithium in the atmosphere. The detection and strength 
of lithium absorption is a clear confirmation of the substellar nature of these objects, as no 



Bodenheimer 


1965; 


Nelson et al. 


1993) 



with the still relatively high effect ive tempera t ure representative for the spectral type L this 



i s a cl ear sign of youth < 1 Gyr (jBasril Il998l ; iMagazzu et al.lll993l ). Moreover, iGizis et al. 



( 12 OOOl ) point out that chromospheric activity is primarily observed in older and hence stellar 
mass L dwarfs. They also find a clear anti-correlation between the presence of lithium and 
H Q emission. Thus, early L dwarfs with a very weak or non-existent H a emission line are 
more likely young brown dwarfs rather than stars. 

The final sample was compiled out of ~ 250 L dwarfs cataloged in D. Kirkpatrick's online 
archival as of 2004 (start of this project) and consists of 12 young L dwarfs which fulfilled 
the following requirements: 



• the object is within 30 pc to the Sun 

• the Li I absorption at 6708 A is present and therefore a clear confirmation of their 
substellar nature, and indication of youth. 

• lack of strong chromospheric activity (H Q emission) which is preferentially observed in 
stellar mass L dwarfs. 

• the objects are isolated in the sense of no known close companion so far. 



These age and proximity criteria helped to enlarge the range of mass and separation over 
which the survey would be sensitive and to maximize the discovery probability. The proper- 



3 http: / /www. DwarfArchives.org 
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ties of all selected targets are summarized in Table [TJ Although there are no Lithium measure- 
ments available for the two mid-L dwarfs 2MASSW 004521+1634 and 2MASSI 0835-0819, 
they were included since they belonged to the 10 closest L dwarfs known and because of 
their brightness compared to other L dwarfs of the same spectral type. A similar decision 
was made for Kelu-1 whi ch had only a weak detection of the Lithium absorption line but was 



apparently overluminos (IMartm et al.lll999t iLeggett et al. 



2001 



Golimowski et al.ll2004f) and 



showed photometric and spectroscopic variability (j Clarke et al. 



2002 



20031 ) . Furthermore, 



the nearby T dwarf 2MASSW 055919-1404 was added to calibrate the spectral differential 
imaging method for HST/NICMOS. 



4. Observations and SDI data reduction 

The observations of the twelve isolated L dwarfs and the calibration T dwarf were 
scheduled as HST program GO 10208 (PI: W. Brandner) from September 2004 to July 2005. 
The sources were observed with the NICMOS1 (NIC1) camera, providing a high-resolution 
pixel scale of 0'.'0432 and a field of view (FoV) of 11" x 11", in the two narrowband filters 
F108N (1.08 /im) and F113N (1.13 /xm). The observation sequences consisted of 4 exposures 
at two different detector positions in each filter. This two - point dither pattern facilitates an 
optimized PSF sampling, effectively rejects bad pixels and gains redundancy against cosmic 
ray events. Altogether, two orbits per target were necessary to achieve the required S/N ratio 
of 6 - 8 for a limiting magnitude of J = 20 mag in the narrowband filters, and a brightness 
difference of ~ 6 mag at a separation of 0.3". All data were acquired in MULTIACCUM 
mode and the total integration times per filter were 2560s (F108N) and 2816s (F113N), 
respectively. In order to limit possible HST breathing variations in the PSF, differential 
imaging observations in the F108N and F113N filters were obtained in each of the orbits. 

For the analysis of the data sets, the HST pipeline reduced images were used, followed by 
an additional bad pixel masking in order to compensate for the non-optimal bad pixel mask 
available in the pipeline. The SDI reduction was then accomplished with a custom made IDL 
program: First, a square aperture around the brown dwarf is extracted from all four images in 
each filter. To make use of the full integration time, the extracted frames were then co-added 
after they were aligned with respect to the first image of the related filter. This alignment 
was done with a 2D cross-correlation and a FFT-shift based combination, which provides 
a sub-pixel accuracy. In the next step, a sub-pixel re-sampling was performed in order to 
transform the observation in the F108N filter to the same A/D scale as in the F113N filter 
(where A is the respective observing wavelength and D is the HST primary mirror diameter). 
The final difference between the two filters was then calculated after additional alignment of 
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the two master images with the shift algorithm, followed by a flux calibration to account for 
any flux losses during the previous re-scale process. In cases where a clear positive residual 
remained in the reduced image, the difference was also obtained separately for the individual 
dither positions (hence only two images per filter). This provides two independent detections 
of a real signal, or the exclusion, if it only exists at one position and is most likely caused 
by a cosmic ray remnant or hot pixel buried in the central core of the PSF. 



5. Results 

The visual examination of the reduced images reveals that for nine of our targets the 
residuals are at the level expected from photon noise and no significant positive signal can be 
detected. The final images of these data sets are shown in Figures [2] and |3j The remaining 
single positive signals in some of the targets in Figure [2] just close to the original PSF center, 
are examples for residuals caused by bad pixels which could not be corrected by the additional 
bad pixel masking, since they were buried in the bright core of the PSF. They are not visible 
in both of the "single" position SDI images simultaneously and thus do not correspond to a 
real signal of a possible faint companion. For the four brightest and closest sample targets 
2MASSW 0045+1634, LSR 0602+3910, 2MASSI 0652+4710 and 2MASSI 0835-0819 (Figure 
[3]) the final SDI results show systematic mirrored dark and bright features in the ~ 6 - 7 
central pixels of the images. These are common residuals in contrast-limited cases and are 
most probably caused by the systematic variations in the PSF structure due to differential 
aberrations of the filters. Since the other five targets are much fainter and their residuals are 
mostly dominated by photon noise, this pattern is no longer visible. If they were brighter 
their residuals would look the same. 

To asses the overall sensitivity for the survey, the standard deviation of the background 
residuals after PSF subtraction was calculated for each data set as a function of radial 
separation from the host brown dwarf. The limits are then based on the median of the 
statistical errors at the various radii and the primary F108N magnitude. Figure H] shows the 
derived 3a narrowband detection limits. The bold solid line represents the mean of the data 
sets, while the shaded area corresponds to the standard deviation. At a separation of 0.05" 
from the primary, a mean contrast of AF108N = 3.5 mag can already be achieved. For the 
brighter brown dwarfs in the survey we reach a contrast as large as ~6.7 mag at separations 
> 0.4" from the primary. For the fainter brown dwarfs in the sample we still reach a contrast 
of ~5.4 mag for separations >0.4". In order to determine the lower mass limit sensitivity 
of the survey, one can translate the contrast (in combination with the brig htness of the 



primary) into a mass using the theoretical COND evolutionary models from iBaraffe et al. 



-9 - 



( 120031 ). First of all, we c onverted the know n J- b and magnitudes of the p rimary L dwarfs 
into Y- band magnitudes. Hillenbrand et al.l ( 120021 ) and lHewett et al.l (120061 ) showed that the 
Y — J colors remain relatively constant at 1.0 ±0.15 mag from early-L through late-T type 
dwarfs. Thereafter, we applied the mass-luminosity relationships of the Y- band (I. Baraffe, 
private communication) which are closest to the F108N bandpass. Figure [5] illustrates the 
resulting mass limits for 1 Gyr (left image) and 0.5 Gyr (right image). Similar to Figure HJ 
the bold solid line corresponds to the mean value of the individual observing runs and the 
shaded areas represent the standard deviation. 

Altogether, the smooth results after the SDI reduction indicate the absence of a plane- 
tary mass companion with a minimum mass of 6-11 Mj up at an age of 1 Gyr and between 
0.07" and ~ 1.4" around the brown dwarfs (~ 1 - 40 AU, depending on their distance to the 
Sun). This is an area usually most strongly affected by PSF residuals for non-differential 
imaging observations. For primaries with an age of 0.5 Gyr our observations even exclude any 
planetary mass object down to 5 - 7Mj up . This result validates our observing and data anal- 
ysis strategy. For the three remaining L dwarfs, the observations revealed very interesting 
results which will be discussed in the following sections. 



5.1. Two newly resolved brown dwarf binaries 

Even though we tried to avoid known binaries during the sample selection process, 
the observations revealed that two of the targets are close brown dwarf binary systems. 
Figure [6] shows the NICMOS detection images where both systems are clearly resolved in 
two components. 



5.1.1. Kelu-lAB 



Kelu-1 was discov ered to be a binary system by iLiu fc Leggettl (120051 ) and also inde- 
pendently resolved by iGelino et al.l ( 120061 ). both with AO from the ground and just four 
months before our own HST observation in July 2005. This solved the by then unexplained 
overluminosity of the "single" object. The HST observation presented here, is now the first 
high-resolution confirmation of the physical association of both components, by proofing 
that they are a common proper motion pair. The system showed clear evidence for orbital 
motion with a significant increase in separation by 15 mas in four months out to 299.8 ± 
0.2 mas and a position angle (PA) of 221.33° ± 0.04°, hinting at a short orbit binary. This 
makes Kelu-1 AB a good target for a monitoring program. The extensive follow up obser- 



vations derived with the Very Large Telescope (VLT) and all results concerning this system 
are discussed in Stumpf etal. (120091 ; 2010b, submitted). 



5.1.2. 2MASSW 0310+1648 AB 



The second discovered very low mass binary is 2MASSW 0310+1648 AB. With an origi- 
nal spectral type of L9, this brown dwarf belongs to the still peculiar L/T transition objects 
which exhibit several unusual characteristics. Resolved for the first time as a binary system 
with almost equally bright components, it adds another system to the already higher binary 
frac tion among the brown dwarfs that span the transition between the L and T spectral class 
(e.g. lBurgasser et al.ll2006l |2010| ). The first astrometric measurements yield a separation of 
204.3 ± 0.4 mas and a PA of 206.4° ± 0.1° on 200 4 Sept ember 24. Including follow up 
observations, this target is discussed in lStumpf et al.l (120101 ) . 



5.2. A planetary companion candidate 

The HST observations of the L4.5 brown dwarf 2MASSW 033703-175807 (hereafter 
2M0337-1758) revealed a clear residual signal in the final SDI reduced image, close to the 
center of the brown dwarf. Unlike for all other L dwarfs in this survey, this signal is also 
visible in the two individually reduced detector position images. Therefore, it cannot be 
caused by a single residual bad pixel. Further, with a S/N ratio of 7 at its peak in the 
final image, it is well above the noise level. Figure [7] displays the resulting images with the 
single detector position results in the middle and the final result, including all observation 
exposures, at the bottom. 



5.2.1. Evidence for binarity 

The F108N image of 2M0337-1758 and its final SDI reduced image were used to calculate 
the magnitude difference between the brown dwarf and the companion candidate, as well 
as their separation and position angle. This was achieved by applying aperture photometry 
with the IRAF phot routine in the apphot package. From the identified centroid positions of 
the photometric apertures a separation of 0.087" ± 0.015" (corresponding to 2.52 ± 0.44 AU 
at a distance of 29 pc) and a position angle of 195.6° ± 5.5° was calculated. 

For the photometry, the preliminary assumption that the flux of the companion only 
contributes to the F108N filter and not the F113N filter was necessary, since the result of the 
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applied SDI technique does not provide any information about the actual flux distribution of 
the companion candidate in each of the filters. This approximation is roughly valid for very 
late T dwarfs and even later spectral type objects (see Figured]), but implies that the derived 
magnitude in F108N for the companion is only a lower limit. For further comparison with 
other observations and theoretical models, the instrumental count rates were converted into 
the Vega magnitude scale after aperture correction and using the most recent photometric 
keyword-value as provided by the STScI webpagd^. Using this method, we derived a F108N 
magnitude of 16.54 ±0.02 mag for 2M 0337-1758 A, a brightness difference of 3.88 ±0.03 
mag between the brown dwarf and its companion candidate and hence a F108N magnitude 
of 20.42 ±0.04 mag for the candidate. With a system distance of 29 pc, this corresponds 
to absolute magnitudes of -vlpiosjv = 14.23 ± 0.02 mag and 18.11 ±0.04 mag for A and B, 
respectively. 

Together with the age estimation of the primary brown dwarf, the comparison of the 
derived values wit h the COND e voluti onary models for brown dwarfs and EGPs (extrasolar 



giant planets) by iBaraffe et al.l ( 120031 ) provides a first approximation of the mass and T c g 
for the companion candidate. Since the F108N wavelength range lies much closer to the 
y-band (Ac = 1.03 /mi) than to the J -band, it was assumed that the objects are equally 
bright in F108N and Y. For the following comparison with the isochrones, an extended 
version of the above mentioned evolutionary models including the Y- band (I. Baraffe, private 
communication) was used. A comparison of the derived mass and effective temperature for 
the primary 2M0337-1758 A with those one would derive with its 2MASS J magnitude, 
proves the assumption of M^ 10 8at = My to be acceptable. The result is given in Table |5J For 
the companion candidate we derive a minimum mass of 10- 15 Mj up with a T e ff between ~ 
600 K and ~ 630 K for an age of 0.5 and 1 Gyr, respectively. To investigate how these results 
would change if the companion magnitudes in F108N and Y are in fact not equally bright, 
we can assume two different scenarios as one can see in the right image of Figure [TJ In the 
two relevant bandpasses the flux of the T dwarfs peaks right in the F108N filter and these 
objects will therefore appear brighter in F108N th an in the F - band . Assuming an offset 



of ~0.5 mag for late T dwarfs, similar to the offset iBiller et al.l (120061 ) derived between the 
H- band and their SDI narrow band, the minimum mass of the companion would decrease 
down to 9- 12 Mj up . In contrast, for a planetary mass object with 5 Mj up at 1 Gyr the water 
absorption band extends into the F108N filter and this object would appear dimmer than 
in the y-band. Here, the offset might not be as large but to derive a valuable upper error 
limit, we use the same value. This results in an increase of the minimum mass to 11-17 Mj up 
for ages of 0.5 and 1 Gyr, respectively. The results for the companion candidate are also 



4 http://www. stsci.edu/hst/nicmos/performancc/photometry/postncs_keywords. html 
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summarized in Table [21 including the larger uncertainties due to the different assumptions. 



As of fall 2009, there are only nine T dwar fs known, which are classified as T8 or 
latejfl and only five of them have a T R ff ~ 600 K (Warren et al. 20071; Delorme et al. 2008 



Burningham et al. 2008 ; Burgasser et al. 2008 ; Burningham et al. 20091 ) . Hence, the com 



panion adds up to a still very small number of just recently detected very cool brown dwarfs 
and the mass estimation places it right at the theoretical mass limit of ~ 13 Mj up for deu- 
terium burning, which is mostly used as the boundary to distinguish between brown dwarfs 
and giant planets. Nevertheless, the dominant uncertainty in translating the derived magni- 
tude difference to brown dwarf and even planetary masses, is the uncertainty of the evolution- 
ary models, which are not yet fully calibrated by observations, especially the mass-luminosity 
relation as a function of age (see further discussions below). 



5.2.2. A possible background object? 

To clarify if the detection might be caused by a background object, the position of 
2M0337-1758 on previous HST/WFPC2 images was analyzed. These images were taken 
on 2000 December 25 in the broadband filter F814W and the medium band filter F1024M. 
No source could be detected in the WFPC2 frames above the noise level (corresponding to 
m F8i4W = 24 mag) at the 2004 NICMOS position of the brown dwarf, but the comparison 
of the two positions allowed the first proper motion determination for 2M 0337-1758 with 
/i a cos<5 = 0.247 ±0.02"/yr and fi s = - 0.024 ±0.028"/yr. This result will also be crucial to 
confirm the common proper motion nature of the binary in follow-up observations. 

Analyzing the probability of a background source from a more physical side, we can 
already exclude a faint, stellar background source (like a M dwarf). These objects have, 
like L dwarfs, a relatively flat continuum between 1.08 - 1.13 fxm and would be therefore not 
detected by our SDI technique. On the contrary, this method does not exclude a distant 
galaxy or a quasar. Based on the fact that we have a point-like remaining signal with an 
approximate magnitude of J ~ 19.5 mag, the candidate is unlikely a "normal" galaxy and 
even too bright to be an active galactic nuclei (AGN). However, we cannot exclude a quasar 
or other peculiar background object (e.g. Wolf- Rayet galaxy). Therefore, the next step will 
be to establish common proper motion of the two components and to exclude any background 
object hypothesis. 



5 http://www. DwarfArchives.org 
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Discussion 



The direct detection of planetary mass objects still remains a challenging endeavor. 
In the past few years many systematic high-contrast direct imaging surveys for planetary 
mass companions around a large number of you ng, nearby stars have been conducted with 



HST ( Luhman et al. 20051; Lowrance et al. 20051) and with AO at large ground -based tele- 



scopes (e.g. 



Masciadri etaJ 120051 iMetchev! [200Q : lLafreniere et all I2OO7I : iBiller et all 12007 



Kasper et al.ll2007al . the latter two employing SDI). However, only very few planetary mass 
or very low mass brown dwarf companions have been discovered in these surveys. The first 
directly imaged planetary mass compani on was detected around the young brown dwarf 
2MASS 1207-3932 (hereafter 2M1207) by Ichauvin et aP (bool l2005ah with an estimated 
mass of 5-8 Mj up and a separation of ~ 55 AU. Shortly after, several more planetary 
mass o bject detections were r eported to yery young primaries with an age < 60 Myr: GO 



Luhman et a 



Lupib flNeuhauser et al.ll2005f) . AB PicB ( Chauvin et al.ll2005bl). D H Tau B (lltoh et al 



200 6b]), Oph 1622-2 405 B flJavawardhana fc Ivanovi 12 006 



CHXR73B flLuhman et all l2006bf ) and 1RXS J1609-2105 B (lLafreniere et al.l 120081 ) . With 



Close et a 



2005 



20071 ) 



mass estimates between ~7-20 Mj up (thus at or just above the brown dwarf/planet bound- 
ary) and projected separations of 55-300 AU from their primary, their nature is still exten- 
sively debated and the exact number of planetary mass companions remains controversial. 
Such large separations suggest a formation process more similar to brown dwarfs (via frag- 
mentation) rather than giant planets (via core accretion, out to a maximum distance of ~ 
10-15 AU; or via gravitational disk instability which demands a very high massive initial 
disk). Alter natively, planets cou ld form close-in to the star and then migrate outward to a 
larger radii. Martin et al. (|2007 ) show t hat such an outward migration could in principle be 
possible. Recently, iMarois et al.l ( 120081 ) announced the first directly imaged multiple plan- 
etary system (HR8799b,c,d) with mas ses between ~ 7- 1 Mj up at separations between 24 



and 68 AU. At almost the same time, iKalas et al.l (12008) det ected the so far lowest mass 
companion Fomalhautb (~ 3 Mj up ) and lLagrange et al.l ( 120091 ) claimed the discovery of the 
closest companion ever directly imaged to a star (/3Pictorisb with ~ 8 AU). While Fomal- 
hautb is again at a relatively large separation, HR8799b,c,d as well as /3Pictorisb, are the 
first directly imaged companions around a star, whose formation could be explained by the 
common core accretion or disk instability scenarios. 

But not only the formation process hinders the exact assignment of all these companions 
to the planetary regime. So far, the masses of these planet candidates are entirely based 
on evolutio nary models. These models have been developed in signific ant detail over the 
past years (IBurrows et al.lll997l ; lBurrowsll2005l ; iBaraffe et al.ll2003l . l2008l ) but remain purely 
theoretical, especially for planets and very young ages due to the lack of direct detections 
of well - characterized giant exoplanets. For example, the original mass estimate for GQ 
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Lupb ranged from 1 to 42 Mj up (INeuhauser et al.ll2005l ). This wide range in mass estimate 
is primarily caused by the use of different theoretical evolutionary models (with different 
initial conditions, see lBurrows et al.lll997l ; lBaraffe et al.ll2002tlWuchterl &: Tscharnuterll2003l ) 
to determine the mass as a function of luminosity, temperature and estimated age. The 
resulting huge discrepancy shows the difficulty of determining an absolute mass for such 
directly imaged low mass objects and especially at such young ages. 

As mentioned in §[TJ the first surveys for brown dwarf binaries already revealed that 
the resolved systems showed a tendency toward a high mass ratio distribution (M 2 /Mi = q 
> 0.8) with a clear peak at q~l. Further high-resolution surveys on ultracool dwarfs 
searching for binary systems strengthened this trend (see Figure[8] on the left). Although 
an underestimation of the number of low mass ratio systems cannot be excluded due to the 
lack of sensitivity for q < 0.5 in most of these surveys, the clear peak to equal mass bin ary 
systems is not an observational b ias effect (fo r a detailed review see iBurgasser et al.ll2007l and 
for an extensive Bayesian study lAllenl 120071 ) . The same surveys unfolded a semi-major axis 
distr ibution which pe aks around 3-5 AU with only ve ry few separations beyond ~ 15 AU 



Bouv et al. 



of Umbreit et al. 



2003 



Table 1 in 



2005 and Allen 



2003). 



Burg asser et al.ll2007l ; supported by numerical simulations 



In the following, the results of our HST survey will be compared with respect to the 
findings discussed above. The resolved and confirmed brown dwarf binary Kelu-1 AB with 
its mass ratio q = 0.82 and a separation of 5.59 AU at the time of discovery, as well as 
2MASSW 0310+1648 AB with its mass ratio q ~ 1 and a separation of 5.17 AU, are in 
full agreement with the previous findings and are an important contribution to the brown 
dwarf binary statistics. The non-detection of any companion in this survey at separations 
larger than 10 AU is also consistent with the separation statistics. The general observed 
binary fraction of very low mass stars (VLMS) and brown dwarfs, as found in direct imaging 
surv eys is ~15-20%, and might be a s high as ~25-30% (for extended discussions see 
e.g. IBurgasser et all 120071 ; lAllenl 120071 ) . With a total sample of 13 brown dwarfs, the two 



confirmed binary systems yield a binary fraction of 15l^ % ( or ^-13 % ^ 2M0337-1758 AB 
is confirmed as a third binary system^) for our survey. This is consistent with the results 
from previous VLMS/BD surveys. 

The discussion about the companion candidate to 2M0337-1758 can only be preliminary 
until second epoch observations confirm the companionship. Yet, it shows some very inter- 



6 The uncertainty of the measured binary fraction was estimated by searching for the fraction of binaries 
at which the binomial distributi on P b (x,n,p) for x positive events in n trials falls off to 1/e of its maximum 
value (cf. iBasri fc Reinersll2006l ). 
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esting properties. With a preliminary mass ratio of 0.18-0.22 (depending on its estimated 
age), this system would add to the very small number of resolved brown dwarf binaries with 
q < 0.5. One of them, which has a resolved and confirmed brown dwarf pri mary and a sim- 



ilar mass ratio is 2M1207AB (q w 0.3 if the most actual mass estimation of iMohanty et al. 



(120071 ) is taken into account). But even if its secondary has a mass just below the deuterium 
burning limit, and therefore appears to be similar to the possible 2M0337-1758 AB system, 
the 2M1207AB system has an unusual large separation and is much younger (~8 Myr). At 
these very young ages the surveys suggest a flatter mass ratio distribution than that for field 
or older cluster binaries (see shaded bins in Figure [8] on the left), thus the low q value might 
not be too unusual in the case of 2M1207AB. Another interesting system in this account 
is SCR 1845 AB. This M8.5-T6 field binary (where the secondary was discov ered during a 



SDI su rvey as well) reveals also a relatively low mass ratio with q « 0.41 (IKasper et al. 



2007bl ). Additionally, the SCR 1845 AB system is, with a separation of ~4.5 AU, similar 
tightly bound as 2M0337-1758 AB. Hence, both systems might be potential candidates for a 
dynamical mass determination on a reasonable time scale. This will further help to calibrate 
evolutionary models under a different aspect. The mass ratio distribution in binary systems 
exhibits a clear correlation with the primary mass (see Figure [8] on the right). While G 
dwarf primaries show a relatively flat mass ratio distribu tion, the ratio for M dwarfs, VLMS 



and brown dwarfs peaks more towards unity (see also iBergfors et al.l |2010| and references 
therein). Although a q value of 0.18-0.22 is, independent of age and separation, still much 
larger than the mass ratios known from general star - giant planet systems (determined with 
RV, transit etc. methods), it does not completely exclude a brown dwarf- PMO system. Due 
to the lower initial mass of the primary the mass ratio of such a system would be a priori 
higher than for a main sequence star -giant planet system. Thus, such a q value for our 
system probes an interesting ratio regime which thus far has not been explored for brown 
dwarfs. 

The minimum mass estimate of 10-15 Mj up for the companion candidate of 2M0337- 
1758 does not a priori exclude a classification as a planet. Ne v erthe less, a comparison of 



four of the ~ 600 K T dwarfs (all field dwarfs) in iLeggett et al.l (120091 ) shows, that at least 
three of them (so far single objects) have similar mass estimates, just right below or above 13 
Mj up . Therefore, 2M 0337-1758 is more likely a brown dwarf - brown dwarf binary system 
with a so far unusual very low mass secondary which possibly straddles the T/Y transition, 
rather than a system including an exoplanet. Follow-up observations of the system are 
proposed to achieve a definite classification. Even more important, a larger number of 
further detections of objects in this mass and temperature regime are necessary to reliably 
calibrate the theoretical models. 
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7. Conclusions 

While previous direct imaging surveys focused especially on searches for planetary mass 
objects around young stars in nearby associations, we obtained high-resolution observations 
with HST/NICMOS of 12 L dwarfs within 30 pc of the sun, using the spectral differential 
imaging technique in the two narrow-band filters F108N and F113N. To this date, this is the 
only SDI survey for planetary mass companions around field brown dwarfs (not related to 
any young association or moving group) conducted from space or ground-based telescopes. 
Further, it is the only one which uses the spectroscopic feature of water absorption instead 
of the methane feature. 

The survey resolved two brown dwarf binaries (Kelu-1 AB and 2MASS 0310+1648AB) 
for which common prope r motion has already been confirmed in follow-up observations 



flStumpf et all I2009I . |2010[ ). The overall binary fraction of 15+^ %, as well as separations 
of the binaries smaller than 6 AU are consistent with previous brown dwarf binary statistics. 
The mass ratios of q > 0.8 confirm the preference for equal mass systems as found in previous 
direct imaging surveys. 

Furthermore, tentative evidence was found for a very low mass companion close to the 
planetary mass regime around 2MASS 0337-1758. It has a very low T eff « 600-630 K and 
an unusual mass ratio (g ~ 0.2) compared to the vast majority of previously found brown 
dwarf binaries. Follow-up observations are proposed to confirm common proper motion of the 
system and to derive more physical properties like colors for spectral type classification. In 
addition, these new observations will provide second epoch astrometry for the determination 
of orbital parameters and hence an initial dynamical mass estimation. If gravitationally 
bound and depending on the age, the 2MASSW 0337-1758 system will be one of the closest 
resolved brown dwarf binaries (0.087"± 0.015", corresponding to 2.52 ± 0.44 AU at a distance 
of 29 pc) with the so far coldest companion ever directly imaged and the least massive 
companion to any L or T dwarf. Therefore, including its multiplicity status, it will then be 
an important testbed in the newly explored < 700 K temperature regime and might imply 
new constraints on the existing formation scenarios. 
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Fig. 1. — Observed spectra of L and T dwarfs (from S. Leggett's online archive) with the superimposed 
bandpass of the HST F108N and F113N filters. The T dwarf plot a dditionally includes the synthetic 
spectrum of a 5 Mj up , 1 Gyr old planetary mass object with T e ff = 312 K (jBurrows et al.1120031 ). For a better 
comparison all flux is arbitrarily scaled. 
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2MASSWJ0051 1 0-1 544 rebinned & oversampled 




2MASSWJ01 0332+1 935 rebinned & oversampled 




2MASSIJ08251 9+21 1 5 rebinned & oversampled 




Fig. 2. — Result images for five of the nine L dwarfs without any companion: The image on the left 
represents the added image of all 4 exposures taken in the F108N filter with the upper cut-level reduced for 
a better visibility of the background (thus not representing the real peak values). The right image displays 
the final reduced image after the subtraction of the two filters. The positive signals visible in the central 
part of the final images of 2MASS 0051-1544 and 2MASS 0825+2115, are not simultaneously present in each 
of the individual reduced detector position images and are therefore residuals from bad pixels buried in the 
original PSF. The overall smooth result images confirm the correctness of our reduction procedure. 



- 25 - 




-26 - 



2MASSWJ004521 +1 634 



30 
25 
20 
15 
10 
5 




30 
25 
20 
15 
10 
5 




0.6 
0.5 

0.4 ! 
i 

0.3 J 
; 

o.2 ; 

O.i 

o.o 



rebinned & oversampled 



5 10 15 20 25 30 
pixel 

LSR0602+3910 



ft 



0.8 

0.6 ! 
i 

i 

0.4 1 

i 
( 

0.2 



5 10 15 20 25 30 
pixel 

2MASSIJ0652+4710 




5 10 15 20 25 30 
pixel 

2MASSIJ0835-0819 





30 






25 






20 


_ _ 


pixel 


15 






10 






5 












0.6 
0.5 

0.4 8 
t/i 

0.3» 
c 

0.2° 
O.i 



5 10 15 20 25 30 
pixel 





30 




25 




20 


03 




X 
Q_ 


15 




10 




5 








0.2 
0.1 

u 
B 

00 5 
c 

-0.1 | 
-0.2 



5 10 15 20 25 30 
pixel 



rebinned & oversampled 




5 10 15 20 25 30 

pixel 



rebinned & oversampled 




5 10 15 20 25 30 
pixel 



rebinned & oversampled 




0.10 



0.05 £ 



0.00 



-0.05 



10 15 20 25 30 
pixel 



Fig. 3. — same as Figure [2] These are the four brightest and closest targets of the sample. The systematic 
mirrored dark and bright features are most probably caused by the systematic variations in the PSF structure 
due to differential aberrations of the filters. However, no other significant positive signal can be found in the 
vicinity of the brown dwarfs. 
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Fig. 4. — The image illustrates the achievable sensitivity in magnitude difference AF108N for our SDI 
survey as a function of angular separation. The bold solid line corresponds to the mean 3tr detection limit 
of all observing runs (excluding the 2 confirmed and clearly resolved binary systems). The upper and lower 
shaded curves correspond to the standard deviation. 
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Fig. 5. — Comparison of achieved detection limits for companion masses at ages of 1 Gyr (left image) and 
0.5 Gyr (right image). The deri ved contrast limits o f each observing run were converted into masses based 



on the evolutionary models from lBaraffe et aL ( 20031 ) as described in the text. The images display the mean 



detectable mass limit with the shaded areas representing the standard deviation. 




Fig. 6. — NICMOS images of the two resolved L dwarf binaries: on the left the bright, initially classified 
L2 dwarf Kelu-lAB, and on the right the much fainter L9 dwarf 2MASS 0310+1648 AB 
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Fig. 7. — The reduced images for the L4.5 dwarf 2MASSW 0337-1758. The top image represents the added 
image of all 4 exposures in filter F108N and the bottom one the residual image of all exposures after the 
SDI reduction (same as in Figures [2] and [3]). In the middle panel, the left inlay displays the residual image 
at position 1 on the detector, while the right displays the residuals at position 2. The clearly remaining 
positive signal is visible in all three residual images and thus any contamination by a bad pixel can be ruled 
out. The brightness difference in the final image is AF108N = 3.88 mag. The orientation of all inlays is the 
same and the north direction is indicated by the arrow. 
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Fig. 8. — The left image illustrates the mass ratio distribution computed for the currently known VLM 
binary systems, based on the online VLM Binaries Archive (www.vlmbinaries.org) as of July 2009. Spec- 
troscopic binaries without known mass ratios were excluded. The distribution shows clear evidence of a peak 
near unity. The shaded bins represent the systems with age estimates < 10 Myr, which have on average 
larger separations compared to the older binary systems. The right image displays the compari son of mass 
ratio distributions for G stars, early M stars and ultracool dwarfs, as adopted from I Alien! (|2007l) . The mass 
ratio appears to increas with decreasing mass of the primary. 



Table 1. List of selected Targets 



Object name 


RA 


DEC 
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Distance a 
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2MASSW 033703-1758 


03 h 37 m 03.6 s 


-17°58'07.9" 


L4.5 


29 


8 


23.6 


15.6 


2004 Sep 24 


LSR 0602+3910 


06 h 02 m 30.5 s 


+39°10'59.2" 


LI 


10.6 


7 


20.8 


12.3 


2005 Apr 15 


2MASSI 065230+4710 


06 h 52 m 30.7 s 


+47°10'34.8" 


L4.5 


11.1 


14 


21.4 


13.5 


2005 Mar 24 


2MASSI 082519+2115 


08 h 25 m 19.6 s 


+21°15'52.0" 


L7.5 


10.7* 


10 


22.5 


15.1 


2005 Feb 12 


2MASSI 083542-0819 


08 h 35 m 42.6 s 


-08°19'23.7" 


L5 


8.3 


? 


21.0 


13.2 


2005 Jun 14 


2MASSW 124646+4027 


12 h 46 m 46.8 s 


+40°27'15.1" 


LI 


25.1 


11 


22.3 


15.0 


2004 Oct 22 


Kclu-1 


13 h 05 m 40.2 s 


-25°41'06.0" 


L2 


18.7* 


1.7 


22.1 


13.4 


2005 Jul 31 


2MASSW 152322+3014 


15 h 23 m 22.6 s 


+30° 14' 56.2" 


L8 


18.6* 


9 


24.4 


16.3 


2004 Sep 30 


2MASSW 055919-1404 


05 h 59 m 19.1 s 


-14° 04' 48.9" 


T4.5 
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Ref. 



3,4 
1 
1 
1 
1 
1 
5 

1,2 
5 
1 

2,6 
1 



CO 



• denotes lacking observations 

no observation since atomic lithium forms into LiCl for T c ff < 1400 K l|Burrows fc Shard Il999l ') and objects that cold are 
brown dwarfs by default 

a Spectrophotometric distance estimates from the given references unless otherwise noted. Distance error is ~ 10% 

* Distance from trigonometric parallax in the given references 
b from the CDS Simbad service 

c from 2MASS All-Sky Point Source Cataloge 

d Uncertainties of apparent magnitudes are ~ 0.1 mag 



References. 
f5) ICruz et ail ||2003 



Dahn et alj \20<d±. (3 l lWilson et all l|2003l ,. (4') ISalim et all (|2003h . 
Burgasser et all i2003al) 
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Table 2. Mass estimations for the companion candidate of 2MASSW 0337-1758 



Object 


My [ mag ] 


0.5 


Gyr 


1 


Gyr 






M Jup 


Teff [K] 


M Jup 


Teff [K] 


Brown Dwarf" 


14.23 ±0.02 


54 ±2 


1870 ±60 


68 ±1 


1945 ±35 


companion candidate 6 


18.11 ±0.50 


> io 


600 ± 70 


15 ±2 


630 ± 70 



Note. — Masses and effective temperatures derived from: 
a DUSTY00 evolutionary models bv lchabrier et all fepQoh : 
b COND03 evolutionary models by iBaraffe et aD J2003h 
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